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CO : ABSTRACT 

a\ : 

The electroweak phase transition in a constant hypermagnetic field is stud- 
■ ied in the Standard Model. The symmetry behaviour is investigated within the 

consistent effective potential of the scalar and magnetic fields at finite tempera- 
^ | ture. It includes the one-loop and ring diagram contributions. All fundamental 

fermions and bosons are taken into consideration with their actual masses. The 
only free parameter is the Higgs boson mass which is chosen to be in the en- 
Q\ • ergy interval 75 GeV < m# < 115 GeV. It is found that for the field strengths 

H ~ 10 22 — 10 23 G the electroweak phase transition is of first order but a baryo- 
genesis condition is not satisfied. For stronger fields it becomes of second order. 
Hence it is concluded that the smooth hypermagnetic field does not generate the 
strong first order phase transition and the baryogenesis does not survive in the 
,iZ< ■ Standard Model. The comparison with the results of other approaches is done. 

> 

•rH , 

^ ' 1 Introduction 

d; 

Among nowadays problems of high energy physics there are two ones which, 
at first glance, are not connected with each other. These are the value of the 
Higgs boson mass m# and the magnetic field strengths H which can be present 
in the early universe (see surveys ||T|, [Q]). They are of paramount importance 
for particle physics and cosmology. For instance, a large scale homogeneous 
hypercharge magnetic field could essentially influence the type of the electroweak 
(EW) phase transition making it strong first order 0, ||. An interest to this 
problem increased recently when it has been realized that without external fields 
a standard baryogenesis does not hold in the minimal Standard Model (SM). 

In Refs. m, ||, 0, @ the influence of the constant hypermagnetic field on 
the EW phase transition has been investigated. In the former one the EP was 
computed in a tree approximation and the result that the presence of Hy makes 
the weak first-order phase transition stronger has been derived. In Ref. |j| the 
temperature dependent part of the EP was calculated in one-loop order whereas 
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the external field has also been allowed for in tree approximation. Therein it was 
found that for the field strengths Hy > 0.3 — 0.5T 2 , where T is temperature at 
the transition, the standard baryogenesis survives. As this paper is concerned, 
we would like to notice that for the weak first-order phase transition the fluctu- 
ations are essential, the one-loop approximation to the EP is not trusty and the 



correlation corrections must be included |J, jTIJ, fll| . Moreover, no investigation 
of the EP curve in strong fields has been carried out in Ref. ||. In fact, it was 
just assumed that the phase transition is of first order and the jump of the order 
parameter is as at the zero external field. Furthermore, the role of fermions at 
high temperature and strong fields has not been investigated. The same has rel- 
evance to the papers ||, ||, fTj. But as it will be seen in what follows, fermions 
significantly influence a vacuum dynamics in the environment. 

To make a link between recent studying of symmetry behaviour in the ex- 
ternal hypermagnetic field and the already obtained results for the case of usual 



magnetic field fll2 |, |15| , we notice that in the broken phase Hy is connected 



with H by the relation H = Hy cos where 6 is the Weinberg angle. So, all 
investigations dealing with symmetry behaviour in the magnetic fields at high 
temperature are relevant to the considered case of Hy in the respect of the form 
of the EP curve at different T,Hy. The hypercharge field influences the scalar 
field condensate at tree level, as magnetic field in the Higgs model, and acts to 
restor symmetry. That was the reason why it has been allowed for in lower order 
in Ref. 0. But, as we will see for strong fields and heavy m#, the form of the EP 
curve in the broken phase is very sensitive to the change of the parameters. In 
the restored phase, there is a number of terms having order ~ (gHy) 3 ^ 2 T which 
can influence the temperature of the phase transition. So, to have an adequate 
picture of the phenomenon investigated the radiation corrections in the field must 
be calculated in both the broken and the restored phases. 

In the present paper the EW phase transition in the constant strong hy- 
percharge magnetic field is investigated within the consistent EP including the 
one-loop and ring diagram contributions. All bosons and fermions are taken into 
account with their actual masses (in particular, the t-quark mass is 175 GeV). 
So, the only free parameter remains the mass m#. We assume it to be in the 
energy range 75 GeV < m# < 115 GeV, in order to take account of the mod- 
ern experemental low limit m# > 90 GeV. We calculate the contribution of ring 
diagrams in the external fields. As it is well known, these diagrams cancel the 
imaginary terms of the one-loop effective potential and the total potential is real 
at sufficiently high temperatures |27| . Due to the ring diagrams in the field, there 
is also cancellation of an instability generated in strong magnetic fields in the W- 
boson sector. This instability appears because of the presence in the VF-boson 
spectrum of the tachyonic mode e 2 = p| + M 2 — eH (see survey |j35"| ). This mode 
is the transversal one. So, to treat the problem carefully the ring diagrams of 
the mode with the transversal effective mass at nonzero H, T have to be added. 
This requires calculation of the W- boson mass operator at high temperature and 
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strong fields. We present the relevant results also. With such the term included 
the total EP is real at sufficiently high temperatures and suitable to investigate 
symmetry behaviour. 

For the bosonic part of the Salam- Weinberg model the phase transitions in 
magnetic fields at high temperature have been studied in one-loop approximation 
in Refs. 0- @, |5f, jz|, p[. However, the aspect of the EW phase tran- 



sition, which was not investigated, is the influence of the correlation corrections 
described by the ring diagrams at high temperatures and strong fields. At zero 
field it was studied in detail in Refs. |27||, f9j, [|10|] but for not heavy mass m t . In 



particular, in Ref. m t was chosen of order ~ 110 GeV. So, for the present day 
experimental data it should be revised. 

In what follows, considering EP in the broken phase we will write H for the 
usual external magnetic field remembering that it equals to H = HycosO. The 
Z-component of the field Hy is screened by the scalar field condensate <p c . The 
constant extenal field is a good approximation for the description of the inital 
stage of the first order EW phase transition when the bubles are not large, as it 
was discussed in Ref. [[|. 

The content is as follows. In Sects. 2, 3 the one- loop contributions of bosons 
and fermions to the EP V^ 1 ) (T, H, C ) are calculated in the form convenient for 
numerical investigations. In Sect. 4 we compute the contributions of ring di- 
agrams. In Sect. 5 the EP for the restored phase is calculated. In Sect. 6 
symmetry behaviour at high temperatures and strong external fields is inves- 
tigated and it is shown the EW phase transition is of first order for the field 
strengths H ~ 10 22 — 10 23 G but the baryogenesis condition is not satisfied. For 
stronger fields it becomes of second order. Thus, we come to conclusion that 
in the SM baryogenesis does not survive under smooth external hypermagnetic 
field. The comparison of our results with that of other approaches and discussion 
are given in Sects. 6, 7. 

2 Boson contributions to V^(T, H, (j) c ) 

The Lagrangian of the boson sector of the Salam- Weinberg model is 



4 a 4 



L = — A F« V F<? - -jG^GT + (D^)+(D^) 



(1) 



where the standard notations are introduced 

_abc a b 



F« u = d,A a u -d u A«+ge abc A^A c u , 
G^ u = d il B v — duB^, 

D, = d, + UgA«r a + -ig'B,. (2) 
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The vacuum expectation value of the field $ is 

The fields corresponding to the W-, Z-bosons and photons, respectively, are 

W± = ^=(Al±iAl), 

Z» = 7 ^=,( 9 Al- 9 >B,), 

A ^=^p(9 >A l + 9B,). (4) 

To incorporate an interaction with an external hypermagnetic field we add the 
term ^HHy to the Lagrangian. The value of the macroscopic magnetic field 
generated inside the system will be determined by minimization of free energy. 
Interaction with a classical electromagnetic field is introduced as usually by split- 
ting the potential in two parts: = + A R , where A R describes a radiation 
field and A = (0, 0, Hx 1 , 0) corresponds to the constant magnetic field directed 



along the third axis. We make use of the gauge- fixing conditions [35 



d,W^ ± ieApW** =F i 9 ^^ = C^x), (5) 

- ^ + gf 2 ) 1 ^. = C M , (6) 

where e = gsinO, tangO = g'/g, 4> , 4> z are the Goldstone fields, ^, are the gauge 
fixing parameters, C ± , C z are arbitrary functions and C is the value of the scalar 
condensate. In what follows, all calculations will be done in the general relativistic 
renormalizable gauge (||), (||) and after that we set = choosing the unitary 
gauge. 

To compute the EP in the background magnetic field let us introduce 
the proper time (s-representation) for the Green functions 

oo 

G ab = -i f dsexp{-isG- iab ) (7) 
b 



and use the method of Ref. [28], allowing in a natural way to incorporate 



the temperature into this formalism. A basic formula of Ref. |28[ connecting 
the Matsubara-Green functions with the Green functions at zero temperature is 
needed, 

+oo 

Gf {x, x'- T) = J2(-l) {n+[x])ak Gt b (x - [x]Pu, x' - npu), (8) 
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where Gf 1 is the corresponding function at T = 0,/? = l/T,u = (0,0,0,1), the 
symbol [x] means the integer part of x±j /3, = 1 in the case of physical fermions 
and <Tfc = for the boson and the ghost fields. The Green functions in the right- 
hand side of formula (§) are the matrix elements of the operators Gk computed 
in the states | x', a) at T = 0, and in the left-hand side the operators are averaged 
in the states with T ^ 0. The corresponding functional spaces U° and U T are 
different but in the limit of T — ► U T transforms into U°. 

The one-loop contribution into EP is given by the expression 

V {1) = -^Tr log G ab , (9) 

where G ab stands for the propagators of all the quantum fields W , (/> ,-■■ in the 
background magnetic field H. In the s-representation the calculation of the trace 



can be done in accordance with the formula [24] 

ds 



f (1 s 

Tr log G ab = - —trexpi-isG^ 1 ). (10) 
o 

Details of calculations based on the s-representation and the formula (H) can be 



found, for instance, in Refs. | 28fl , [12]] , p8[ . The terms with n = in Eqs. (^|), (|9|) 
give the zero temperature expressions for the Green functions and the effective 
potential V^ 1 ), respectively. They are the only terms possessing divergences. 
To eliminate them and uniquely fix the potential we make use the following 
renormalization conditions at H, T = | |12| : 

(11) 



d 2 V(<p,H) 1 

T^p | H =0,0=5(0)- 2> 

dV&H) 



H=0,cp=8(0)= 0, (12) 



d 2 V^,H) , . 2 



\H=0,<p=6(0) z 



m 1 |, (13) 



where V((p, H) = + + • • • is the expansion with respect to the num- 
ber of loops and 5(0) is the vacuum value of the scalar field determined in tree 
approximation. 

It is convenient for what follows to introduce dimensionless quantities: h = 
H/H (H = Ml/e) A = <p c /S(0),K = m 2 H /Ml B = (3M w ,t = 1/B = 
T/M W ,V = Vj Hi and M w = §5(0). 

After the reparametrization the scalar field potential is explicitely expressed 
in terms of the ratio K, 

V(°) = ^ + ^m^(-0 2 + ^). (14) 
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Remind that h is the electromagnetic component of the hypercharge field hy 
which is unscreened in the broken phase. In the restored phase it will be conve- 
nient to work in terms of the initial fields and we will carry out the corresponding 
calculations later. 

The renormalized one-loop EP is given by the sum of the functions 

V, = V (0) + V (1) (0, h, K) + cu {1) (0, h, K, r), (15) 

where is the one-loop EP at T = 0, which has been studied already in Ref. 



It has the form: 

V (1) =V£ + V« (16) 



where 



- ±-hHog(2h)} + ^[-2A 2 + (ft 2 + h4?)log(h + tf) 

+ (h 2 - h4> 2 )log \h-4> 2 \}+ i^ah((f) 2 - h)d(h - 4> 2 ), (17) 



v£>= k Si n 2 e w (-<P 2 + l<P 4 ) 

+ 5^(1 + 2_ 3 4 2) 

4vr v 2cos 4 # A 2 r yr 4 r r; 

aK 2 r/ 9 , 4 , 2 1,, . 30 2 - 1 . 27 , 4 21 , 2l , . 

and cj« is the temperature dependent contribution to the EP determined by the 
corresponding terms of formulae (|8|), (Q) with 

We outline the used calculation procedure considering the W-boson contribu- 
tion as an example |38|, 



2tt 



ll^/^ii^f^^v^. (19) 

1 

As Eq. fllTD ) this expression contains an imaginary part for h > <f) 2 appearing due 
to the tachyonic mode £ 2 = p 2 + M 2 — eH in the V^-boson spectrum [R3], [CT, 



3(J. It can be explicitly calculated by means of an analytic continuation taking 
into account the shift s — > s — iO in the s-plane. Fixing <p 2 /h > 1 one can rotate 
clockwise the integration contour in the s-plane and direct it along the negative 
imaginary axis. Then, using the identity 



and integrating over s in accordance with the standard formula 

b , 



/ dss* 1 - 1 exp( as) = 2(-) n/2 K n (2v^6), (21) 

J s a 



(a, b > 0), one can represent the expression (|i"9|) in the form 



i?eo;« = -4-4(3^ + u) X - u 2 ), (22) 

7T D 



where 



oo oo 



^ = E E -^(nftrp); x p = (0 2 + h + 2p/0 x / 2 , (23) 



p=0n=l n 



OO 



uj. = £ V -Ki{nBy), y = {<f - h) 1 ' 2 . (24) 
We have in the range of parameters <fi 2 < h after analytic continuation 



oo 



uJ 1 = --J2 LJ ^Y 1 (nB\y\), (25) 



00 2 



^ 2 = T -K 1 (nBz),z = (0 2 + h) x l\ (26) 

K n (x), Y n (x) are the Bessel functions. The imaginary part of u$ is given by the 
expression 

h 00 I v I 

J/nwi = -2a- £ ^ Ji(nS | y |), (27) 

n=l n 

Ji(x) is the Bessel function. As it is well known, the imaginary part of EP is 
signaling the instability of the system. In what follows we shall consider mainly 
symmetry behaviour described by the real part of the EP. As the imaginary part, 
it will be cancelled in the consistent calculation including the one-loop and ring 
diagram contributions to the EP. 



Carrying out similar calculations for the Z- and Higgs bosons, we obtain [[12 

U * = _6_ E 2/3 2U2 K ^ fl)> ( 28 ) 

7T cos^ v w n z B z cos v 

tee t 2 



« E ^^(nS 1 t |) 

n=l 
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where the variable t = [K w (^^-)] 1/2 at 3(f) 2 > 1 and series with the function 
Y 2 (x) has to be calculated at 30 2 < 1. The corresponding imaginary term is also 
cancelled as it will be shown below. 

The above expressions (|16"D , (0), (|28|) , (|29|) will be used in the numerical 
studying of symmetry behaviour at different H, T. There is a cancellation of a 
number of terms from the zero-temperature contributions given Eqs. ([16]) and 
T-depended ones. This fact has a general character and was used in checking of 
the correctness of calculations. 



3 Fermion contributions to V^\H,T,(/) C ) 

To find the convenient form of the fermion contribution to the EP let us consider 



the standard unrenormalized expression written in the s -representation |33 

V p =" (-l) n / ^e- {m2 f s+p2n2/As \eHs)coth{eHs), (30) 

° n n=-oo ^ S 

ruf is a fermion mass. Here, we have incorporated the equation (||) to introduce 
the temperature dependence. In what follows, we shall take into account the 
contributions of all fermions - leptons and quarks - with the values of masses 
equalled to the present day data. Usually, considering symmetry behaviour with- 
out external fields one restricts themself by the t-quark contribution, only. But in 
the case of the external field applied this is not a good idea, since the dependence 
of on H is the complicate function of the parameters mj, eH, T. At some 
fixed values of H, T fermions with the definite corresponding masses are domi- 
nant. For instance, at high temperature the liding term of Vf is ~ H 2 log^j. 
Hence it follows that light fermions are important . In general, a very complicate 
dependence on the field takes place. We include this in the total, carrying out nu- 
merical calculations and summing up over all the fermions. Now, separating the 

+oo oo 

zero temperature terms by means of the relation = 1 + and introducing 

— oo 1 

the parameter Kf = m 2 /M 2 = G\ ukawa i 1 g' ! , we obtain for the zero temperature 
fermion contribution to the dimensionless EP, 

V/(M) = £E^(- 2 ^ 2 + ^ 4 -^W) 

a 2 h 2 2 I q f I h. 

rE(9? ¥ log -^ 



TT 



a E^iogr^-p^-) + (2C'(-1) - \)q)h 



2 



7T y 1 2 | qf | h 6' 

+ I^V + ti^-i^U/l^log^F, (31) 



S 



where qj is a fermion electric charge, the sum = J2 (leptons) + 3 J2 (quarks) 

f f=i f=i 
counts the contributions of leptons and quarks with their electric charges. The 

function Ti is defined as follows |3l| (see also survey ||35| ): 

logr x (a:) = fdylogr(y) + ~x(x-l)-~xlog(2Tr). (32) 
o 



The finite temperature part can be calculated in a way described in the pre- 
vious section. In the dimensionless variables it looks as follows: 

+ ((2P + 2)h + K^Bn)] } (33) 

art J > 

Again, a number of terms from Eqs. ( |3"TD and ( |3"BD are cancelled being summed 
up, as in the bosonic sector. 

These two expressions and the boson contributions obtained in Sect. 2 will 
be used in numerical investigations of symmetry behaviour. 



4 Contribution of ring diagrams 

It was shown by Carrington |J that at T ^ the consistent calculation of the 
EP based on generalized propagators, which include the polarization operator 
insertions, requires that ring diagrams have to be added simultaneously with the 
one- loop part. These diagrams essentially affect the phase transition at high 



temperature and zero field [p7[|, H, [110]. Their importance at T and H ^ was 



also pointed out in literature IfLql , [jig] but, as far as we know, this part of the 
EP has not been calculated, yet. 

As it is known PU , [g7|, the sum of rine diagrams describes a dominant con- 
tribution of large distances. It non-negligibly differs from zero only in the case 
when massless states appear in a system. So, this type of diagrams has to be 
calculated when a symmetry restoration is investigated. To find the correction 
V r i n g(H,T) at high temperature and magnetic field the polarization operators of 
the Higgs particle, photon and Z-boson at the considered background have to 
be computed. Just these calculations have been announced in Refs. |15| , [ p~6f| . 
Then, V ring (H,T) is given by series depicted in figures 1, 2. Here, a dashed line 
describes the Higgs particles, the wavy lines represent photons and Z-bosons, the 
blobs represent the polarization operators in the limit of zero momenta. As it is 
also known M, in order to calculate the contribution of ring diagrams not the 
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Figure 1: The Higgs field ring diagrams giving contribution to the effective po- 
tential. Blobs stand for the neutral scalar field polarization operator calculated 
at zero momentum. 




Figure 2: The photon and Z-boson ring diagrams giving contribution to the effec- 
tive potential. Blobs stand for the polarization operators of the fields calculated 
at zero momenta. 



total polarization operators n^fc, T, if) but only their limits at zero momenta, 
noo(^ — 0,T, if), are sufficient. This limit, called the Debye mass, can be cal- 
culated from the EP of the special type. This fact considerably simplifies our 
task. 

Now, let us turn to calculations of V ring (H,T). It is given by the standard 
expression P7| , ||, JT5[ : 



Vring = "^^{[^(^ + n oo(0)] 3/2 - M 3 (0)}, (34) 

where trace means the summation over all the contributing states, M(0) is 
the tree mass of the corresponding state. The functions n 00 (0) are: n 00 (0) 
=H(k = 0, T, if) for the Higgs particle; n 00 (0) = n 00 (&; = 0, T, if) - the zero-zero 
components of the polarization functions of gauge fields in the magnetic field 
taken at zero momenta. The above contributions are of order ~ g 3 (A 3 / 2 ) in the 
coupling constants whereas the two-loop terms are to be of order ~ g 4 (X 2 ). For 
n o(0) the high temperature limits of polarization functions have to be substi- 
tuted which have order ~ T 2 for leading terms and ~ g<fi c T, (gH) l l 2 T((l) c /T « 
1, (gHf/ 2 /T « 1) for subleading ones. 

For the next step of calculation, we remind that the effective potential is 
the generating functional of the one-particle irreducible Green functions at zero 
external momenta. So, to have 11(0) we can just calculate the second derivative 
with respect to <p of the potential V^'(H, T, </>) in the limit of high temperature, 
T >> (p,T >> (eif) 1//2 , and then set = 0. This limit can be obtained by 
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means of the Mellin transformation technique (see, for instance, ||38|| ) and the 
result looks as follows: 



C 



f A2 



6 



+ 



an 



2cos 2 6„ 



+ 



+ 



[^(3A0 2 -5 2 (O))T 2 



g ^3 T _ «, 3A0g-£ 2 (O) 3/2 
cos 3 6> 3 2 



2vr 



-e#T( 



eif) 



1/2 



(35) 



3 J2 Gf determines the fermion contribution of 



3 

The parameter Cf = J2 G\ 

i=l i=l 

order ~ T 2 having relevance to our problem. We also have omitted ~ T 4 contri- 
butions to the EP. The terms of the type ~ log[T/ /(</>, H)] cancel the logarithmic 
terms in the temperature independent parts (|15|), (|30|). Considering the high tem- 
perature limit we restrict ourselves to linear and quadratic in T terms, only. 

The one else important expression, which also should be taken into account, 
is the linear in H term of the zero temperature EP Eq. (0), which looks as 
follows: 



qjH 



(36) 



It significantly influences symmetry behaviour and contributes to the Debye mass 
in strong fields. 

Now, differentiating these expressions twice with respect to <p an d setting 
4> = 0, we obtain 



n,(o) 



&vv>{<i>,h,t) 



1 0=0 



6e 2 



+ 



+ 



24/3 
2a x ^ 



+ 



3e 2 



sin 2 26„, sin 2 6. 



\- 3/3 2 ' 
a/2 



q f H | K f 



Sir sin 6 w p 2 2 



(37) 



The terms ~ T 2 in Eq. (57) give standard contributions to temperature mass 
squared coming from the boson and the fermion sectors. The if-dependent term 
is negative since the difference in the brackets is 3v^2~C(~ 



I I' 

2' 2< 



- 1 



-0,39. 



Formally, this may result in the negativeness of the 11(0)0 for very strong fields 
(eif) 1 / 2 > T. But this happens in the range of parameters where asymptotic 



11 



axpansion is not applicable. Substituting expression fl3"7| ) into Eq. (|34] ) we obtain 
(in the dimensionless variables) 

As one can see, the last term of this expression cancels the fourth term in the Eq. 



(|35|), which becomes imaginary at 3(f) 2 < 1. This is the important cancellation 
preventing the infrared instability at high temperature. 

Before we proceed, let us note that Eq. (^) contains the other term (the last 
one) which becomes imaginary for strong magnetic fields or small 2 . It reflects 
the known instability in the H^-boson spectrum which is discussed for many years 



in literature (see papers |}I2"|| , |i~5|j , fll6] , p8| and references therein). But it also 
will be cancelled out when the contribution of ring diagrams with the unstable 
mode is added. 

To find the if-dependent Debye masses of photons and Z-bosons the following 
procedure will be used. We calculate the one-loop EP of the W-bosons and 
fermions in a magnetic field and some " chemical potential" , fi, which plays the 
role of an auxiliary parameter. Then, by differentiating them twice with respect 
to fi and setting /i = the mass squared m 2 D will be obtained. Let us first 
demonstrate that in more detail for the case of fermion contributions where the 
result is known. 

The temperature dependent part of the one-loop EP of constant magnetic 
field at a non-zero chemical potential induced by an electron-positron vacuum 



polarization is p3 |: 



V /e 1 L = ^E(- 1 ) m j / -^exp(^--m 2 s)(eHs)coth(eHs)cosh((3l^ (39) 

1=1 o 

where m is the electron mass, e = gsinO is the electric charge and the proper-time 
representation is used. Its second derivative with respect to \i taken at ji = can 
be written in the form, 

= ™0>£- / -exp(-m 2 s - (3 2 l 2 /As)coth(eHs). (40) 

Expanding coth(eHs) in series and integrating over s in accordance with formula 
(pTh we obtain in the limit of T » m,T » (eH) 1 / 2 : 



B-ljW^iWI) + l(&M Kl {m0l) + •••]. (41) 

The series in / can easily be calculated by means of the Mellin transformation (see 
Refs. , JT^[). To have the Debye mass squared it is necessary to differentiate 
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Eq. (HOI) with respect to 8 and to take into account the relation of the parameter 



/i with the zero component of the electromagnetic potential: /i — > ieA [y|. In 
this way we obtain finally, 

m 2 D = g 2 sm 2 d\^- - ^m 2 + 0((m/3) 2 , {eH8 2 ))] . (42) 

This is the well known result calculated from the photon polarization operator 
[29[|. As one can see, the dependence on H appears in the order ~ T~ 2 . To find 



the total fermion contribution to m 2 D one should sum up the expression ( S2|) over 
all fermions and substitute their electric charges. 

To find m 2 D for Z-bosons it is sufficient to allow for the fermion coupling to 
the Z-field. It can be done by substituting u — > i(g/2cos9 + gsin 2 9) and the 
result differs from Eq. (|42| ) by the coefficient at the brackets in the right-hand 
side which has to be replaced, g 2 sin 2 9 — > g 2 ( 4 + tang 2 6). One also should add 
the terms coming due to the neutral currents and the part of fermion-Z-boson 
interaction which is not reproduced by the above substitution: 

m ^ = sik (1+4sm ^ )T2 - (43) 



Now, let us apply this procedure to the case of the Pl^-boson contribution. 
The corresponding EP (temperature dependent part) calculated at non-zero T, u 
in the unitary gauge looks as follows, 



K? } = - G ^- 2 f:]- 2 ex P {-m 2 s-l 2 p 2 /As) (44) 

l=i n 



3 

+ 4sinh(eH s)]cosh(/3l{i) . 



sinh(eHs) 



All the notations are obvious. The first term in the squared brackets gives the 
triple contribution of the charged scalar field and the second one is due to the 
interaction with a VT-boson magnetic moment. Again, after differentiation twice 
with respect to /i and setting \i = it can be written as 

d 2 V^ eH o2 d ™°fds . m 2 s l 2 p 2 eH.. 3 A . fA . 

- w n = — -8 2 t—; V / — exp( — — A ) tt -+Asinh(s) . (45) 

dy? | M=0 2n 2H df3 2 ^J o s yK eH 4s n sinh(s) v n v 1 

Expanding sinh~ l s in a series over Bernoulli's polynomials, 

i — s co r> 

E§(-2^) fe , (46) 



sinks s r-~L k\ 

fe=0 
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and carrying out all the calculations described above, we obtain for the W- boson 
contribution to m 2 D of the electromagnetic field, 

m 2 D = 3g 2 sin 2 9[^T 2 - ^T(m 2 + gsin6H) 1/2 - — {gsin9H) 

+ 0(m 2 /T 2 ,(gsin6H/T 2 ) 2 )]. (47) 

Hence it follows that spin does not affect the Debye mass in leading order. Other 
interesting feature is that the next-to-leading terms are negative. 

The contribution of the H^-boson sector to the Z-boson mass m 2 D is given by 
the expression ( |47|) with the replacement g 2 sin 2 9 — > g 2 cos 2 8. 

Substituting the expressions (^) and fl4"T|) into Eq. (0), we obtain the photon 
part V^g, where it is necessary to express masses in terms of the vacuum value 
of the scalar condensate C . In the same way the ring diagrams of Z-bosons 
V r z ing can be calculated. The only difference is the mass term of Z- field and the 
additional term in the Debye mass due to the neutral current ~ v^i^vZ^. These 
three fields - <p, 7, Z- which become massless in the restored phase, contribute 
into V r i n g(H, T) in the presence of the magnetic field. At zero field, there are also 
terms due to the H^-boson loops in Figs. 1, 2 . But when H 7^ the charged 
particles acquire masses ~ eH and can be neglected. 

In the restored phase, the IV-bosons do not interact with the hypermagnetic 
field and therefore give no field dependent contributions. 

A separate consideration should be spared to the tachyonic (unstable) mode in 
the H^-boson spectrum: p\ = p\ + M 2 — eH. First of all, we notice that this mode 
is excited due to a spin interaction and it does not influence the G 00 (k) component 
of the jy-boson propagator. Secondly, in the fields eH ~ M 2 the mode becomes 
a long range state. Therefore, it should be included in V ring (H,T) side by side 
with the other considered neutral fields. But in this case it is impossible to take 
advantage of formula ( |34"|) and one has to return to the initial EP containing the 
generalized propagators. 

For our purpose it will be convenient to make use of the generalized EP written 
as the sum over the modes in the external magnetic field [ITB], [TEfl: 



+00 



TT +OO „ J OO 

^S = ^E/^ E MW + i^ + n^tf))], (48) 

Z7r P l=-oo_ oc Z7F n=0,<r=0,±l 

where to t = ^ , e 2 n = pj + M 2 + (2n + 1 - 2a)eH and U(H, T) is the radiation 
mass squared of jy-bosons in a magnetic field at finite temperature. Denoting as 
Dq 1 (p 3 , H.T) the sum tuf + e 2 , one can rewrite eq. ([48]) as follows: 
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+ S £ / ^09[l + ("? +P3 + — eH)-'U(H, T)) 

P l=- 00 _ 00 

+ J2 log[l + D (e 2 n ,H,T)U(H,T)}}. (49) 

Here, the first term is just the one-loop contribution of VT-bosons, the second one 
gives the sum of ring diagrams of the unstable mode (as it can easily be verified 
by expanding the logarithm into a series). The last term describes the sum of 
the short range modes in the magnetic field and should be omitted. 

Thus, to get V™ng table one has to compute the second term in Eq. ( f49|) . In the 
high temperature limit we obtain: 

pH 

V r ZT ble = ^{( M - - eH + m,T)) 1 ' 2 - (M 2 W - eH) 1 / 2 }. (50) 

By summing up the one-loop EP and all the terms V r i ng , we arrive at the total 
consistent in leading order EP. 

Let us mention the most important features of the above expression. It is seen 



that the last term in Eq. (50) exactly cancels the "dangerous" term in Eq. 
So, the EP is real and no instabilities appear at sufficently high temperatures 
when II(if , T) > — eH. To make a quantitative estimate of the range of 
validity of the total EP it is necessary to calculate the H^-boson mass operator 
in a magnetic field at finite temperature and hence to find U(H,T). This is a 
separate and enough cogent problem which is considered in detail in a separate 
publication. Here, we only adduce the result of H(H, T) calculations p4|| : 

U unstable (H, T) = < n = 0, a = +1 | \ n = 0, a = +1 > 

= a[2Q,9Q(eH) 1/2 T + iA(eH) 1/2 T], (51) 

where the average value of the mass operator in the ground state of the M^-boson 
spectrum | n = 0, o = +1 > was taken. This expression has been computed in 
the limit eH/T 2 << 1,B = M W (H,T)/T << 1, which is a good approximation 
since, as it will be shown below, typical inverse temperatures for the symmetry 
restoration are B ~ 0.1 — 0.3. Side by side with the real part responsible for the 
radiation mass squared the expression (|5lD contains the imaginary one describing 
the decay of the state. Its value is small as compare to the real part and of order of 
the usual damping constant at high temperature. So, Im H(H, T) can be ignored 
in our problem. The radiation mass squared is positive and acts to stabilize the 
spectrum. At H = no screening is produced in one-loop order, as it should be 
at finite temperatures for transversal modes |17||- Thus, we see that at high tem- 
peratures the effective H^-boson mass squared, (M^,) e jj = M 2 — eH + U(H,T), 
is positive. Therefore, no conditions for H^-boson condensation discussed in Refs. 

P5| are realized. With this result obtained we conclude that our EP is real 
for temperatures corresponding to the phase transition epoch.*) [] 

1 Expression (|5l|) disagrees with the corresponding one of Ref. [Sol where the average value 
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5 Effective potential of the restored phase 



Having obtained the EP at 7^ we are able to investigate the form of its 
curve in the broken phase and determine the type of the EW phase transition 
for different m#, h. To describe more precise the restored phase one has also to 
calculate radiation corrections to the external hypermagnetic field Hy at high 
temperature. Before doing that let us remind that at = the field Hy is 
completely unscreened. This means that in a covariant derivative describing 
interaction with the external field one should include the U(l)y term: = 

+ \g'B e *\ We set the potential as before, BJf = (0, 0, H Y , 0). 

In the restored phase W-bosons do not interact with H Y . The field dependent 
part of the EP V(4> = 0, Hy, T) is non-zero due to the contributions of fermions 
and scalars. However, the fermion part depends logarithmically on temperature 
(~ (g/ 2 ) JJylogT /T ) and can be neglected as compared to the tree level term 
\H Y . The scalar field contribution to the one- loop EP is 

V£\H Y ,T)= - Ul^iln(TIT ) 

+ MW!l +0(1/r) . (52 ) 

The term logarithmically dependent on T can again be neglected but the linear 
in T part must be retained. Since "hyperphotons" are massless in the restored 
phase we also include the contribution of the corresponding ring diagrams: 

V r Zle d (Hy,T) = -^-[ 2 -(g'/2) 2 T 2 + m 2 Df -MBM!^.-^ /2 )H Y \^ 

} ~V , " (53) 

where m 2 D . = ^g' 2 T 2 Yf is the sum over the fermion contributions to the 

f{R ' L) 

Debye mass of the "hyperphotons" , Yf are the hypercharges of R— and L— lep- 
tons and quarks. Both these expressions have been calculated in a way described 
in the previous sections. 

For convenience of numerical investigations let us rewrite Eqs. ( |52D and (|53D in 
terms of the dimensionless variables h, B: V(H Y , T) Testore d = (H ) 2 v restorec i(h, B), 



^restoredi^h, B*J 



1 h 2 a h 3 / 2 

2cos 2 6 + 3y/2cos 3 6 B 

1 a.7 Ana h 1 ^ 2 h 3 / 2 

3WQcos 2 6B 2 ~ 2V2irBcos6 ~ lQn 2 cosW ' ^ ' 



of the gluon polarization operator in an abelian chromomagnetic field was calculated in a weak 
field approximation and H(H,T) has been found to be zero. Most probably, the discrepancy 
is the concequence of the calculation procedure adopted by these authors when the gluon 
polarization operator was computed at zero external field and then its average value has been 
calculated in the state | n = 0, a = +1 > . Our expression is the high temperature limit of the 
mass operator which takes account of the external field exactly. 
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where a = e 2 /4ir and hy = h/cos9. 



6 Symmetry behaviour in a hypermagnetic field 

Now, let us investigate the EW phase transition in the hypermagnetic field for 
different values of m^. It can be done by considering the Gibbs free energy in 
the broken, G broken (H ext ,(j) c ,T), and the restored, G restored (H Y xt ,T), phases f|, 
pi. The first order phase transition can be determined from two equations: 

GrestorediH^ ,T,0) = Gbroken(H eXt , T, (j)(H eXt ) c ) , (55) 

describing the advantage of the broken phase creation, where 4>{H) C is a scalar 
field vacuum expectation value at given H, T which has to be found as the mini- 
mum position of the total EP, 

dV(H,T,^ = Q 

d<f) c 

Hence the critical field strength can be calculated. In this expression (and below) 
we write for brevity H instead H ext . 

Having obtained the EP in the restored phase, the one-loop EP described by 
formulae (|T6|) , (0), fl28|)-(|33|) and the ring diagram contributions V ring we are 



going to investigate the symmetry behaviour. First, we consider the total EP as 
the function of 2 at various fixed H, T, K and determine the form of the EP 
curves in the broken phase. In this way it will be possible to select the range of 
the parameters when the first order phase transition is realized. After that the 
temperature T c at given field strength (Hy) c will be estimated using Eqs. (|55|). 
©• 

As usually ||35|| , to investigate symmetry behaviour we consider the difference 
V = Re[V(h, 0, K, B) - V(h, = 0, B)] which gives information about symmetry 
restoration. 

In what follows it will be also convenient to express the conditions of the phase 
transition in terms of the dimensionless variables h, B, 0, taking into account the 
relation hy = h/cosO. Then, the Gibbs free energy 

G broken (h ext , 0, B) = y + v'{h, 0, B) - hh ex \ (57) 
has to be expressed through h ext by using the equation 

h^ = h+ dv ' (h g ^ B \ (58) 

where v describes the one-loop and the ring diagram contributions to the EP. 
The phase transiton happens when the relation holds, 

h 2 

—tan 2 6 = v' restored (h, B c ) - v' broken {h, C , B c ). (59) 
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The function v' restored is given by Eq. 
h ext by h. 

The results on the phase transition determined by numerical investigation of 
the total EP are summarized in Table 1. 



54). We also have substituted the field 



h 


K 


T c (GeV) 


<f> c (h,T c ) 




R 


M w (h,B c ) 


0.01 


0.85 


106.47 


0.301662 


0.091 


0.69699 


0.327235 


0.01 


1.25 


122.21 


0.181659 


0.033 


0.36567 


0.230086 


0.01 


2 


145.56 


0.094868 


0.009 


0.16033 


0.186168 


0.1 


0.85 


108.58 


0.275681 


0.076 


0.62459 


0.245186 


0.1 


1.25 


123.54 


0.130384 


0.017 


0.25963 


0.112721 


0.1 


2 


148.39 


0.031623 


0.001 


0.05242 


0.126315 


0.5 


0.85 


108.89 


0.248998 


0.062 


0.56253 


0.49938 i 


0.5 


1.25 


second 


order 


phase 


transition 




0.5 


2 


second 


order 


phase 


transition 





Table 1. 

In the first column we show the hypermagnetic field strength in the broken 
phase (in dimensionless units). In the second and third the mass parameter 
K = mjj/M^ and the critical temperature of the first order phase transition 
are adduced. Next two columns give the local minimum positions <p c (H,T c ) and 
their squard values at the transition temperatures. The last two columns fix the 
ratio R = 246 GeV <p c (h,T c )/T c , determining the advantage of baryogenesis, and 
the 1^-boson effective mass calculated in the local minimum of the EP at the 
corresponding field strengths and the transition temperatures. 

As it is seen, the increase in h makes the phase transition weaker (not stronger 
as it was expected in Refs. [§J, by analogy to superconductivity in the external 
magnetic field). The ratio R is less than unit for all the field strengths, wherease 
the baryogenesis condition is R> 1.2 — 1.5 [19|. Thus, we come to the conclusion 
that external hypermagnetic fields do not make the EW phase transition strong 
enough to produce baryogenesis. Moreover, for strong fields the phase transi- 
tion is of second order for all the values of K considered. These are the main 
observations of our numerical investigations. 

Let us continue the analysis of data in the Table 1. For the field strengths 
H > 0.1— 0.5Hq(Hq = M^/e) the phase transition is of second or weak first-order. 
The W-boson effective mass squared (in dimensionless units) M^,(0 C , h, B c ) = 
4>c(h, B c ) — h + U(h, B c ) is positive for h = 0.01 and h = 0.1. Therefore, the local 
minimum is the stable state at the first order phase transition. For stronger fields, 
when the second order phase transition happens, the effective P^-boson mass 
becomes imaginary. This reflects the known instability in the external magnetic 
field which exhibits itself even when the radiation mass of the tachyonic mode 
is included. But it does not matter for the problem of searching for the strong 
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first order phase transition in the external hypermagnetic field investigated in 
the present paper. The instability has to result in the condensation of W- and 
Z-boson fields at high temperature. 

In Refs. |3 , @ , || , Q it has been determined that the strong hypermagnetic 
field increases the strength of the first order phase transition and in this case 
baryogenesis survives in the SM. Our results are in obvious contradiction with 
this conclusion. To explane the origin of the discrepancies let us first consider 
Refs. [f|, |J where a perturbative method of computations has been applied. 
These authors, considering the EW phase transition, have allowed for the influ- 
ence of the external field at tree level, only. That corresponds to the usual case 
of superconductors in the external magnetic field, and, as a consequence, they 
observed the strong first order phase transition. In fact, the type of the phase 
transition was just ussumed, since no investigations of the EP curve with all 
the particles included for different Hy,T have been carried out. In the former 
paper, the qualitative estimate of the phenomenom considered was given, where- 
ase in the latter one the quantitative analysis in one-loop approximation for the 
temperature dependent part of the EP has been done. Actually, in both these 
investigations the influence of the external field was reduced to consideration of 
the condition (|59|) fixing the transition temperature. The role of fermions and 
VT-bosons in the field was not investigated at all. However, as we have observed, 
the fermions (heavy and light) are of paramount importance in the phase tran- 
sition dynamics. Just due to them the EW phase transition becomes of second 
order in strong fields (for the values of K when it is of first order in weak fields) . 

In Refs. ||, |7| the phase transition was investigated by the method combining 
the perturbation theory and the lattice simulations. As the first step in this 
approach the static modes are maintained in the high temperature Lagrangian. 
The fermions are nonstatic modes and decoupled. So, no reflections of the fermion 
properties in the external fields and, hence, no information on the EP curve could 
be derived in this way. The only fermion remainder is the t-quark mass entering 
the effective universal theory |TD| , [§], 0. In our analysis, it has been observed 
that not only heavy but also light fermions are important in strong external fields 
at high temperature. In fact, for various field strengths the fermoins with different 
masses are dominant and we have taken account of all of them. Moreover, we 
have allowed for all the ring correlation corrections in the external field that also 
influences symmetry behaviour. 

We would like to notice that our perturbative results for the values of K ~ 
0.8 — 0.9 are reliable. They are in agreement with nonperturbative analysis at 
zero field. The external field is taken into account exactly. For these masses of the 
Higgs particles we observed the change of the first order phase transition to the 
second order one with increase in the field strength. The same behaviour takes 
place for K > 1 when parturbative analysis may be not trusty. But, as we have 
discovered, the general picture of the field effects is only quantitatively changed 
for heavy scalar particles. In this case also the first order phase transition in 
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weak fields becomes of second order one for strong fields. These circumstances 
convince us that the assumption of Ref. [3J that the hypermagnetic fields is able 
to make the weak first-order EW phase transition strong enough is not aproved 
by the detailed calculations. 

7 Discussion 

In Refs. || the influence of strong external hypermagnetic field on the EW 
phase transition has been taken into account in tree approximation. Further 
studying of the phenomenon, naturally, has to allow for the radiation and corre- 
lation corrections. This is the problem that we have addressed to in the present 
paper. The main idea was to determine the form of the EP curve in the broken 
phase and find the range of the parameters Hy, K when the EW phase transi- 
tion is of first order. To elaborate that the consistent EP including the one-loop 
and ring diagrams of all the fundamental particles has been constructed. As we 
have seen, the role of fermions and ring diagrams in the external field is crucial 
when the structure of the broken phase is described. The external field was taken 
into consideration exactly through Green's functions. The minimum of the EP 
was found to be stable at sufficently high temperatures when the first order phase 
transition happens. This important property is fulfilled when the ring diagrams of 
the tachyonic mode are included. As a result, no conditions for W- and Z-boson 
condensates are realized at high temperatures at the first order phase transition. 
The condensates could be generated for stronger fields at the second order phase 
transition. But in this case baryogenesis does not survive. 

The influence of strong magnetic fields on the vacuum at high temperature is 
a complicate corporative effect described by the total EP. At some chosen values 
of H,T,K the different terms of it are dominant. To better understand the role 
of fermions in symmetry behaviour let us adduce two terms of the asymptotic 
expansion of the EP in the limit of T — > oo, H — > oo. The first one is the 
term ~ H 2 log^j. Due to this term the light fermions are dominant at high 
temperature. The second term can be derived from the expansion of the zero 
temperature part Eq. (|31|). This expression side by side with the leading term 
~ H 2 log^, which due to a "dimension parameter trading" is replaced by the 

above written term, contains the subleading one ~ —eHm 2 log^- (for details see 

Ref. p2p. This term acts to make "heavier" the Higgs particles in the field. 
As a result, the second order temperature phase transition is stimulated due to 
strong fields. As it is well known at zero field, the correlation corrections relax the 
strength of the first order phase transition |9| - |TT|. This property is provided by 
the structure of the V ring term of the EP (|34]), independently of the field presence. 
The field-dependent terms of the Debye masses of the scalar ([57D , photon (pE7j) 
and Z particles are negative that decreases the mass values. Therefore, the field 
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acts to make weaker the effect of correlations as the zero field case is compared. 
But nevertheless the relaxation effect as such holds. We also have investigated the 
influence of different parts of the EP on symmetry behaviour. It was discovered 
that the change of the phase transition kind with increase in H is due to the 
fermion temperature part of the EP. These remarks help us to have a notion 
about the role of fermions and correlations in strong fields. 

In papers ||, @ the EW phase transition in the hypermagnetic field has 
been investigated by means of the method combining perturbation theory and 
lattice simulations and the main conclusions of Refs. f§], [[| were supported. 
In the former two papers, because of peculiarities of the calculation procedure 
adopted, the effects of the external field due to fermions as well as the correlation 
corrections have not been allowed for. So, from the point of view of the present 
analysis these results also do not reproduce correctly the behaviour of the EP 
curve in the broken phase. 

The values of the Higgs boson mass investigated in the present paper corre- 
spond to the cases when perturbative results are reliable (K = 0.85) and may 
be not trusty (K = 1.25 , 2). However, since the external field is taken into 
account exactly its effects do not depend on the specific K values. As we have 
seen, an increase in Hy makes the EW phase transition of second order for the 
field strengths Hy ~ 0.5 • 10 24 G for all the mass values investigated. For weaker 
fields the phase transition is of first order but the ratio R = 4> C (H,T C )/T C is less 
then unit. Hence we conclude that baryogenesis does not survive in the minimal 
Standard Model in the smooth external hypermagnetic fields. 

The authors thank M. Bordag and A. Linde for interesting discussions and 
remarks and A. Batrachenko for checking of the numeric calculations. One of us 
(V.S.) grateful DFG grant No 436 UKR 17/24/98 for financial support and col- 
leagues from Institute of Theoretical Physics at Leipzig university for hospitality. 
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